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1 Introduction 



In a recent paper, hereafter referred to as [1], we explored the phenomenology of a supersymmetric 
SU{4:)c X SU{2)l X SU{2)r (4-2-2) model [2] with third family Yukawa unification [3, 4]. We showed, 
in particular, that Yukawa unification and neutralino dark matter are readily compatible in this 
4-2-2 model, which should be contrasted with the conclusion reached in [5] within the corresponding 
SO(IO) framework where it is extremely hard to do so. The reason for this difference largely stems 
from the assumption made in 4-2-2 of tree-level gaugino non-universality, which is not possible in 
SO(IO). Thus, in contrast to the latter, there is an additional soft parameter in 4-2-2 which enables 
one, via gluino coannihilation [6], to make third family Yukawa unification and neutralino dark 
matter mutually consistent. In [1] we showed some benchmark points highlighting the characteristic 
sparticle and Higgs mass spectra in the presence of Yukawa unification. A relatively light gluino is 
a characteristic feature of this 4-2-2 [1] as well as SO (10) [5] model. 

In this paper we carry out a more thorough investigation of the 4-2-2 model without insisting 
on Yukawa unification. We assume low scale gravity mediated supersymmetry breaking [7], and 
allow the asymptotic gaugino masses to be non-universal at tree level, compatible with the 4-2-2 
symmetry [8, 9]. We identify a variety of coannihilation scenarios in this more general model which, 
as in [1], also has discrete left-right symmetry [2, 10] (more precisely C-parity) [11]. As in [1], 
non-universal gaugino masses, compatible with the 4-2-2 gauge symmetry, play an essential role in 
realizing gluino [6] and bino-wino [12] coannihilation regions that are consistent with the WMAP 
dark matter constraints. We also explore regions of the parameter space in which the little hierarchy 
problem is partially resolved. We present several benchmark points and the associated sparticle and 
Higgs mass spectra. Many of the new particles are in a mass range which is accessible at the LHC. 

Supplementing 4-2-2 with a discrete left-right (LR) symmetry reduces the number of independent 
gauge couplings in 4-2-2 from three to two. This is because LR symmetry imposes SU{2)i and 
SU{2)r gauge coupling unification condition at Mgut (fi'L = S'r)- In 4-2-2 the matter fields are 
unified into three generations of ip (4, 2, 1), and the antimatter fields are in three generations of ipc 
(4,1,2). If the minimal supersymmetric standard model (MSSM) electroweak doublets come from 
the bi-doublet H(l, 2, 2), the third family Yukawa coupling H tpci^ yields the following relation valid 
at Mgut, namely 

= = Yt = ^Dirac- (1) 

In order to get the correct fermion masses and mixings one possibility is to assume that the SM 
Higgs partially comes from A (15, 2, 2) dimensional representation of 4-2-2 [13]. In other words we 
have to assume that the SM Higgs doublet is a superposition of (1, 2, 2) and (15, 2, 2) Higgs fields. 
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In this case, the superpotential for the Yukawa sector of the 4-2-2 model is 

Yi'iP.ijH + Y^'iP.^pA (2) 

With the MSSM Higgs doublets now arising from a linear combination of the Higgs doublets in 
(1,2,2) and (15,2,2), the Yukawa unification condition in Eq.(l) is in general lost. 

We will assume that due to C-parity the soft mass^ terms, induced at Mqut through gravity 
mediated supersymmetry breaking are equal in magnitude for the scalar squarks and sleptons of the 
three families. The tree level asymptotic MSSM gaugino masses, on the other hand, can be non- 
universal from the following consideration. From C-parity, we can expect that the gaugino masses 
at Mgut associated with SU{2)l and SU{2)f{ are the same (M2 = = Mj"). However, the 
asymptotic SU{4:)c and consequently SU{3)c gaugino masses can be different. With the hypercharge 
generator in 4-2-2 given hjY = ^J^(B — L) + ^J^Isr, where B — L and I^r are the diagonal generators 
of SU{4:)c and SU{2)r, we have the following asymptotic relation between the three MSSM gaugino 
masses: 

Ml = + ^Ma. (3) 

5 5 

The supersymmetric 4-2-2 model with C-parity thus has two independent parameters (M2 and 
M3) in the gaugino sector. The fundamental parameters of the 4-2-2 model are as follows: 

mo, mH^,mHa, M2, M3, Aq, tan/3, sign /x. (4) 

Thus, compared to the NUHM2 model of [14], we have one additional parameter in 4-2-2 which 
plays a crucial role in opening up several distinct coannihilation channels each associated with a 
characteristic sparticle and Higgs mass spectrum. 

The outline for the rest of the paper is as follows. In Section 2 we summarize the scanning pro- 
cedure and the experimental constraints that we have employed. In Section 3 we present the results 
from our scan, highlight some of the predictions of the 4-2-2 model, and display some benchmark 
points. Our conclusions are summarized in Section 4. 



2 Phenomenological constraints and scanning procedure 

We employ the ISAJET 7.78 package [15] to perform random scans over the parameter space listed 
in Eq.(4). In this package, the weak scale values of gauge and third generation Yukawa couplings are 
evolved to Mqut via the MSSM renormalization group equations (RGEs) in the DR regularization 
scheme, where Mqut is defined to be the scale at which gi = g2. We do not enforce the unification 
condition = Qi = g2 at Mqut, since a few percent deviation from unification can be assigned to 
unknown GUT-scale threshold corrections [16]. At Mgut? the boundary conditions are imposed and 
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all the (soft supersymmetry breaking) SSB parameters, along with the gauge and Yukawa couplings, 
are evolved back to the weak scale Mi- The impact of the neutrino Dirac Yukawa coupling in the 
running of the RGEs is significant only for relatively large values (~ 2 or so) [17]. In our model we 
expect the largest Dirac coupling to be comparable, at best, to the top Yukawa coupling (~ 0.5). 

In the evaluation of Yukawa couplings the SUSY threshold corrections [18] are taken into account 
at the common scale MgusY = ^m^^mf^. The entire parameter set is iteratively run between Mz and 
^GUT using the full 2-loop RGEs until a stable solution is obtained. To better account for leading- 
log corrections, one-loop step-beta functions are adopted for gauge and Yukawa couplings, and the 
SSB parameters mj are extracted from RGEs at multiple scales = mj(mj). The RG 15- improved 
1-loop effective potential is minimized at an optimized scale Msusy? which effectively accounts for the 
leading 2-loop corrections. Full 1-loop radiative corrections are incorporated for all sparticle masses. 

The requirement of radiative electroweak symmetry breaking (REWSB) [19] puts an important 
theoretical constraint on the parameter space. Another important constraint comes from limits on 
the cosmological abundance of stable charged particles [20]. This excludes regions in the parameter 
space where charged SUSY particles, such as fi or ti, become the lightest supersymmetric particle 
(LSP). We accept only those solutions for which one of the neutralinos is the LSP and saturates the 
WMAP dark matter relic abundance bound. 

We have performed random scans for the following parameter range: 



< 


mo 


< 10 TeV, 


< 


Ma 


< 1 TeV, 


< 




< 1 TeV, 


-30 TeV < 




<o. 


< 




< 10 TeV 


< 


rriHd 


< 10 TeV, 


1 < 


tan (3 


< 65, 



with /i > , and nit = 172.6 GeV [21]. A more recent estimate by the CDF/DO collaboration 
quotes a slightly larger value rrit = 173.1 [22]. This does not change our conclusions in any significant 
way. 

We first collected 1 million points for the 4-2-2 model. To this we added the results of constrained 
MSSM (CMSSM) and CMSSM with non-universal Higgs (150,000 points each), as these are subsets 
of the 4-2-2 model. All of these points satisfy the requirement of REWSB with the neutralino 
being the LSP in each case. Furthermore, all of these points satisfy the constraint QcDuh"^ < 10. 
This is done so as to collect more points with a WMAP compatible value of cold dark matter relic 
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abundance. After collecting the data, we use the IsaTools package [23] to implement the following 
phenomenological constraints: 



m^± (chargino mass) 


> 


103.5 GeV 


[20] 


rrif (stau mass) 


> 


86 GeV 


[20] 


(stop mass) 


> 


175 GeV 


[20] 


(sbottom mass) 


> 


222 GeV 


[20] 


rrig (gluino mass) 


> 


220 GeV 


[20] 


nih (lightest Higgs mass) 


> 


114.4 GeV 


[24] 


BR{Bs fj.+fj. ) 
2.85 X 10-^ < BR{b 57) 


< 


5.8 X 10-s 


[25] 


< 


4.24 X 10-^ (2a) 


[26] 






O.llll^i (5a) 


[27] 


3.4 X 10-10 <{g- 2)^ 


< 


55.6 X 10-10 (3a) 


[28] 



We apply the experimental constraints successively on the data that we acquire from ISAJET. 
In our plots we exhibit points that satisfy all of the above constraints as well as those that do not 
satisfy the [g — 2)^ constraint. 

3 Results 

The asymptotic gaugino mass parameter Mi is constrained by Eq.(3) at Mqut which is characteris- 
tically different from the asymptotic relation Mi = M2 = M3 in the case of CMSSM. In the CMSSM, 
after RGE running, it is therefore not possible to have the three gauginos nearly degenerate in mass 
at the low scale. Indeed at 1-loop level, we have the relations M3 ~ 6M1 and M2 ~ 2Mi. In the 4-2-2 
model, if we start with a low value at Mgut of the ratio M3/M2, it is possible to realize M3 ~ M2 
at the TeV scale. This scenario offers a relatively light gluino NLSP (next to LSP), so that gluino 
coannihilation plays an important role in determining the LSP relic density [6]. Similarly, if we start 
with a large value of M3/M2 at Mgut, it is possible to have the wino nearly degenerate with the 
bino, in which case one opens up the bino-wino coannihilation channel [12]. This can be seen in 
Figure 1 where we plot results in the M2 - M3, mo - M2 and mo - M3 planes. Shown in gray are 
points that satisfy REWSB and the requirement that the LSP is a neutralino. The other three colors 
depict points that satisfy the above mentioned bounds except the {g — 2)^ constraint. These include 
the WMAP bounds on dark matter relic density and various collider constraints {BR{Bs — > 
BR{h 57), and (s)particle mass bounds). Shown in orange is the gluino coannihilation channel 
which occurs at low M3 and high M2 values. We show the bino-wino coannihilation channel in red 
which occurs at low M2 and high M3. Shown in blue are all other coannihilation channels which 
have been extensively studied in the CMSSM and CMSSM with non-universal Higgs [29]. Figure 2 
is similar to Figure 1 except that we now show the effect of imposing in addition the constraint from 
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{g — 2)^. The allowed parameter space in this case is appreciably shrunk, which is understandable 
because as we increase mg, the MSSM particles gradually decouple from the theory and we end up 
with the SM result for {g — 2)^. We note, however, the two new coannihilation channels introduced 
by the 4-2-2 model are also present in Figure 2. 

In Figure 3 we show results in the Aq - M3/M2, tan (3 - M3/M2, niHu/mHd - M3/M2 and /i - M3/M2 
planes. The green points satisfy all the constraints, while the blue ones satisfy all constraints except 
the one from [g — 2)^. We do not distinguish the bino-wino and gluino coannihilation channels 
in Figure 3. The Aq - M3/M2 plane shows that the {g — 2)^ constraint can only be satisfied for 
\Ao\ < 4TeV. The tan/? - M3/M2 plane shows a wide range of allowed values for tan/3. If we 
focus on the rriHu/i^Hd - M3/M2 plane, we can identify points on the M3/M2 = 1 line with the 
non-universal Higgs model. The CMSSM corresponds to a single point in the mnu/^Hd - M3/M2 
plane, while the mHu/f^Hd = 1 line displays the effect of varying the ratio M3/M2 in the NUHMl 
parameter space [14], with Mi given by Eq.(3). The results in the fi - M3/M2 plane show that in the 
4-2-2 model we have the possibility of realizing relatively small /i values and, hence, there exists the 
potential for ameliorating the little hierarchy problem [30]. 

We exhibit in Figure 4 the various coannihilation as well as the resonance channels by plotting 
the relevant (NLSP) sparticle mass versus the lightest neutralino. The color coding is the same as 
in Figure 3 except that we now suppress the gray background corresponding to the region allowed 
by REWSB and a neutralino LSP. In the m^o(bino) - m^o(wino) plane we show those solutions for 
which the lightest neutralino is essentially a pure bino and the NLSP is primarily a wino. The unit 
slope line in this plane corresponds to the bino-wino coannihilation region. In the m^o (not bino) 
- m^o (not wino) plane we can recognize the possibility for coannihilation between the two lightest 
neutralinos with the lightest neutralino not a pure bino and the NLSP not primarily a wino. This 
case mostly corresponds to the coannihilation of higgsinos. 

We display the gluino coannihilation region in the m^o - rrig plane. Consistent with all known 
bounds, the gluino can be as light as 250 GeV or so. 

We can recognize the A-funnel region in the 2m^o - plane where neutralino annihilation occurs 
through the CP-odd Higgs resonance channel. In the m^o - plane we have the stop coannihilation 
region with the possibility of a stop as light as ~ 200 GeV. Even though the stop coannihilation 
region is found in the CMSSM, the possibility of a light stop consistent with all constraints including 
the one from {g — 2)^ is a rather exciting feature of the 4-2-2 model. The stau-coannilation region is 
displayed in the mj^o - plane. 

In Figure 5 we show the various coannihilation and resonance channels as a function of M3/M2. 
We have plotted on the y-axis mj/m^o in the case of coannihilations and mj/2m^o in the case of 
resonance annihilations of the neutralino, where mj is the NLSP mass. This means that in each 
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figure, we need to focus close to tlie line y = 1 in order to see coannihilations and resonances. The 
M3/M2 - '^x'il^xi pl^^6 shows, once again, that bino-wino coannihilation occurs for a relatively 
large value of M3/M2 > 2.5. The coannihilations shown in this plane close to M3/M2 ~ 1 are due 
to the higgsinos. The M3/M2 - ^g/^xl pl^^^ shows that for gluino coannihilations to work, one 
needs M3/M2 < 0.1 at Mgut- It is, therefore, not surprising that these channels are not found 
in the CMSSM as they need a large splitting in the gaugino sector. The other four planes have 
coannihilation solutions corresponding to M3 ~ M2 as well as regions in parameter space that have 
a splitting between M2 and M3. The M3/M2 - m^/m^o plane is rather interesting as it shows that 
stop-coannihilation prefers a ratio of M3/M2 smaller than 1. This is understandable because with 
an increase in the gluino mass, the stop will become heavier. It seems from the M3/M2 - m/j/2m-jj^o 
plane that annihilation via the lightest Higgs resonance occurs for M3/M2 < 1. However, this may 
just be due to a lack of statistics rather than a trend. 

Finally we present a few benchmark points in Table 1 highlighting phenomenologically interesting 
features of the 4-2-2 model. Point 1 has been taken from the bino-wino coannihilation region while 
point 2 represents a solution in the gluino coannihilation region of the 4-2-2 parameter space. Gluino 
coannihilations play a crucial role in reconciling the neutralino dark matter relic density with Yukawa 
coupling unification in the 4-2-2 model, as discussed in [1]. In Point 3, with a neutralino of mass 
~ 46 GeV, we demonstrate the possibility of resonant enhancement of neutralino annihilation via the 
Z-boson resonance state. Point 4 highlights another distinguishing feature of the 4-2-2 model. In 
the case of the CMSSM, the left-handed selectron is heavier than the right handed selectron because 
of the larger strength of the SU{2) interaction. In 4-2-2, it is possible to achieve mg-^ < mg-^. This 
happens typically for a small M3/M2 ratio. Point 4 corresponds to the bino-wino coannihilation 
channel, while point 5 represents a solution with a relatively small magnitude of /i, so that the little 
hierarchy problem can be ameliorated. 

4 Conclusions 

We have shown that the Higgs and sparticle spectroscopy arising from the supersymmetric SU{4:)c x 
SU{2)l X SU(2)ji model with non-universal gaugino masses can be quite distinct from the well 
studied constrained minimal supersymmetric standard model. Gluino and bino-wino coannihlation 
scenarios in particular, with their associated spectroscopy, are particularly exciting possbilities which 
will be tested at the LHC. 
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Figure 1: Plots in the M2 - M3, mo - M2 and mo - M3 planes. Gray points are consistent with 
REWSB and LSP. Blue, orange and red points satisfy the WMAP bounds on Xi dark matter 
abundance and various constraints from colliders {BR{Bs fi~^fi~), BR{b — > 57), and (s)particle 
mass bounds). Orange points represent the gluino coannihilation channel, red points depict the 
bino-wino coannihilation channel, while blue points represent all other channels. 
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Figure 2: Plots in the M2 - M3, mo - M2 and mo - M3 planes. Color coding same as in Figure 1, 
except that now the orange, red and blue points also satisfy the {g — 2)^ constraint. 
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Figure 3: Plots in the Aq - M3/M2, tan/3 - M3/M2, mHu/mHd - M3/M2 and /i - M3/M2 planes. 
Gray points are consistent with REWSB and x? LSP. Blue points satisfy the WMAP bounds on x\ 
abundance and various constraints from colliders {BR{Bs fi~^fi~), BR{b — > 57), and (s)particle 
mass bounds). Green points also satisfy the {g — 2)^ constraint. 
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Figure 4: Plots in the m^o(bino) - m^o(wino), m^o(notbino) - m^o (not wino) , m^o - rrig, 2m^o - m^, 
rrij^o - mi and m^o - mf planes. Also shown is the unit slope line in each case. Color coding same as 
in Figure 3. 
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Figure 5: versus m^o/m^o, mh/2m^o, mg/m^o, m^/2m^o, m^/m^o and mf/m-^o. Color 

coding same as in Figure 3. 
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Table 1: Sparticle and Higgs masses (in GeV units), with = 172.6 GeV and /i > 0. Points 1, 2 and 
3 respectively correspond to the bino-wino coannihilation, gluino coannihilation and niz resonance 
channels. Point 4 lists the spectrum for a solution for which the left-handed selectron is lighter than 
the right-handed selectron. Point 5 is an example of a scenario in which the little-hierarchy problem 
is ameliorated. Points 1, 2 and 5 satisfy all the constraints listed at the end of Section 2, while points 
3 and 4 satisfy all constraints except {g — 2)^. 



15 



